;***************************************************************
;discription：计算植被的光合有效辐射
;
;param1: nrows(input) 行
;param2: ncols(input) 列
;param3: fpar(input)  入射光合有效辐射的吸收比例，格式为矩阵
;param4: sol(input)   太阳总辐射量，格式目前为单值，非矩阵
;param5: apar(output) 植被的光合有效辐射
;
;author: wangqiang( wensten )  2015-09-14  chengdu
;***************************************************************
PRO PRO_NPP_APAR, nrows, ncols, fpar, sol, apar=apar
  ;
  FOR yindex=0, nrows-1 DO BEGIN
    FOR xindex=0, ncols-1 DO BEGIN     
      ;获取fpar像元值
      fpar_cellvalue = fpar[yindex,xindex]
      sol_cellvalue=sol[yindex,xindex]
      
      IF (fpar_cellvalue EQ -9999.0)||(sol_cellvalue EQ -9999.0) THEN BEGIN
        ;无数据区
        apar[yindex,xindex] = -9999.0
      ENDIF ELSE BEGIN
        ;0.5表示植被所能利用的太阳有效辐射（0.4-0.7um）占太阳总辐射的比例
        apar[yindex,xindex] = sol_cellvalue * fpar_cellvalue * 0.5
      ENDELSE      
    ENDFOR
  ENDFOR
  
END